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1. Introduction 

O 

O In this paper, we are interested in the location of conjugate points along a geodesic 

^ in the volumorphism group P^(M^) of a compact three-dimensional manifold M^, 

O without boundary. As shown in [?], these are typically pathological, i.e., they can 

O occur in clusters along a geodesic, unlike on finite-dimensional Riemannian manifolds. 

O We give an explicit algorithm for finding them in terms of a certain ordinary differential 

^ equation, derived via the WKB-approximation methods of Lifschitz-Hameiri [?] and 

I— I Priedlander-Vishik [?]. We prove that for a typical geodesic in T>^(M^), there will be 

pathological conjugate point locations filling up closed intervals in M"*"; hence typically 
the zeroes of Jacobi fields on the volumorphism group are dense in intervals. 

Interest in the geometry of the volumorphism group dates back to 1966, when 
d Arnold [?] showed that if the volumorphism group is equipped with an Riemannian 

^ metric, then the geodesies of this metric correspond to Lagrangian solutions of the ideal 

fiuid equations. In particular, the right-translated tangent vector to the geodesic will 
^ be an incompressible velocity field satisfying the Euler equation of ideal fiuid mechan- 

^ ics. Arnold then computed the Riemannian sectional curvature of this metric in some 

special cases, showing that it was often negative, and concluded that most fluid flows 
00 were likely to be unstable in the Lagrangian sense. In addition, he found some sections 

^ of positive curvature, and suggested that the existence of conjugate points would imply 

2 that there were some stable perturbations of a fiuid fiow. 

This geometric approach to equations of continuum mechanics has since been used by 
O many authors to find partial differential equations whose solutions represent geodesies 

K> on an infinite-dimensional Lie group. Well-known equations arising in this way include 

^ the Korteweg-de Vries equation, the Camassa-Holm equation, and the Landau-Lifschitz 

^ equation. See for example [?] for an overview. In some cases, this approach only works 

formally: the geodesic equation resulting may not be a genuine ordinary differential 
equation on a Banach manifold in any reasonable topology. On the volumorphism 
group, however, the geodesic equation is very well-behaved in a Sobolev topology: 
Ebin and Marsden [?] proved in 1970 that if we extend the group of volumor- 
phisms T>i^(M) to the space of Sobolev class volumorphisms V^^{M), there is a 
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C°° Riemannian exponential map defined in some open subset of each tangent space 
TidX>^(M). In addition, the differential of the exponential map can be computed in 
terms of the Jacobi equation in the same way as on a finitc-dimcnsional manifold; 
Misiolek [?] proved in 1992 that the parallel transport and curvature operators are 
smooth in the Sobolev topology, so that the Jacobi equation is also an ordinary dif- 
ferential equation. These results allow one to study ideal fiuids using standard Picard 
iteration techniques rather than PDE estimates. 

The first examples of conjugate points on the diffeomorphism group were found by 
Misiolek, who in [?] gave a simple example of a fiuid fiow on having a conjugate 

point, and in [?] gave a more complicated example of a fiuid fiow on P^(T^) having 
a conjugate point. Arnold and Khesin, in their comprehensive guide to the geometry 
of the diffeomorphism group [?], asked whether one could expect conjugate points to 
accumulate along a geodesic, or whether they were isolated as on a finite-dimensional 
Riemannian manifold. (In general, as shown by Grossmann [?], conjugate points on 
an infinite-dimensional manifold may cluster together along a geodesic. Grossmann 
defined monoconjugatc and epiconjugate points in terms of the failure of the differential 
of the exponential map to be injectivc or surjcctive, respectively, and showed that the 
two types need not coincide in infinite-dimensions. In addition, there may be clustering 
of monoconjugatc or epiconjugate points along geodesies.) 

This question was answered by Ebin, Misiolek, and the author [?]. We showed 
that for any two-dimensional manifold without boundary, the exponential map is a 
nonlinear Fredholm map of index zero. Thus epiconjugate points and monoconjugatc 
points are the same, and furthermore all monoconjugatc points have finite order and are 
isolated along a geodesic. On the other hand, we showed that on the three-dimensional 
manifold x (the solid torus with fiat metric), a certain rotational fiow has an 
epiconjugate point that is not monoconjugatc. The author [?] extended this result 
to show that on any three-dimensional manifold, the first conjugate point along any 
geodesic is pathological: either epiconjugate and not monoconjugatc, or monoconjugatc 
of infinite order. That paper used the index form to locate the first conjugate point, 
and the method was based on finding a simple ordinary differential equation at each 
point of the manifold whose solutions would give an upper bound for the location of 
a conjugate point. However, since it relied on the index form rather than the actual 
Jacobi equation, it could not be used to find the location of any conjugate point other 
than the first. 

In this paper, we approximate the Jacobi equation itself by an ordinary differential 
equation at each point of the manifold, by working with highly-peaked perturbations 
of the fiuid fiow. This technique is well-known as the WKB method, first applied in 
quantum mechanics. It was first applied to the linearized Euler equation by Bayly [?] in 
a special case, then more generally by Lifschitz-Hameiri [?] and Friedlander-Vishik [?] 
in the early 90s. By a result of the author [?], the Jacobi equation can be decomposed 
into the linearized Euler equation and a simple linearized fiow equation; as a result, 
we can apply the WKB method to the Jacobi equation as well. We do this in Section 
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[3] to obtain the ordinary differential equation 
|(A(t,x.)^)+..(x.)x^ 



0, 7(0) =0, (7(t),U=0. 



Here A(t, Xo) = Dri[t, Xo)'^Dri{t, Xo) is the metric deformation of the flow r/ at Xo G M, 
while uJo{xo) = curl Uo(xo) is the initial vorticity at Xo- The vector is a fixed unit 
vector in T^^M. We prove that if for some fixed Xo & M and £ T^^M, there is a 



solution 7(t) of (1.1) satisfying 7(0) = and 7(a) = 0, then there is an epiconjugate 
point in D^(M^) located at t = a. 

The freedom to specify any Xo and any G T^^M implies the existence of many 
epiconjugate points. In Section |4} we consider at a fixed Xo & M the set of all a G 



and € S such that (1.1) has a solution orthogonal to and vanishing at t = and 



t = a. By analyzing the equation (1.1) in detail, we show that this set is a collection of 
two-dimensional surfaces in M"*" x S*^, which can be expressed as the graph of a function 
on some open subset of S"^. From this, we conclude that the epiconjugate point locations 
form closed nontrivial intervals in IR+. As we point out with an example, these intervals 
frequently extend to infinity. 

2. Background 

Our basic object is a compact Riemannian manifold M of dimension n > 2, without 
boundary, whose Riemannian metric generates a volume form /i. (For the most part we 
will assume M is three-dimensional, although this is not necessary.) It is a well-known 
result of Arnold [?] that formally the geodesic equation on the volumorphism group 
Vf^{M), consisting of all C°° diffeomorphisms rj: M ^ M with ri*fi = fi, where /i is 
the Riemannian volume form, is the equation of ideal fluid mechanics in Lagrangian 
coordinates, when ©^(M) is given the Riemannian metric 



(2.2) (m o ?7,t> o 77)^2 = / {u,v)or]jj,= / {u,v) jj.. 

JM JM 

Here the tangent space at G T>^{M) is given by 

Tr,V^{M) = {vo ri \ divv = 0}. 
The differentials of the left- and right-translations are given by 

(2.3) DL^{v) = Dri{v) and DRrjiv) = vor]. 



By equation (2.2), the L metric is right-invariant under the group action, but not 
left-invariant unless 77 is an isometry. 

Although the above works formally, to treat ordinary differential equations we need a 
Hilbert manifold structure, which we obtain by enlarging both T>^{M) and T^X'^(M) to 
the spaces of Sobolev-class diffeomorphisms and if* vector fields, with s > n/2 + 1. 
Ebin-Marsden established that the geodesic equation for a diffeomorphism rj is a smooth 
second-order differential equation on this manifold, so that it has smooth solutions and 
in particular a smooth exponential map exp;^: Q C Tidl^^(M) V^^i^M), defined by 
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expi^(Mo) = 77(1), where rj is the geodesic with 77(0) = id and 77(0) = Uo- (By right- 
invariance of the metric, it is no loss of generahty to assume all geodesies begin at the 
identity map, which we will do from now on.) If M is two-dimensional, then exp;^ is 
actually defined on all of TidP^(M) (a corollary of a famous result due originally to 
Wolibner [?]). If M is three-dimensional, it is still notoriously unknown whether expj^j 
can be defined on the entire tangent space. 

By right-invariance, the geodesic equation may be decomposed into two equations: 
the flow equation and the Euler equation. The flow equation relates material variables 
to spatial variables, and is given by 

(2.4) ^(t,x) =M(t,r/(t,x)). 
The Euler equation for the velocity fleld u is 

(2.5) ^ + P{VuU) = 0, 

where P is the orthogonal projection of an arbitrary vector fleld onto the space of 
divergence-free vector flelds, given explicitly by 

(2.6) P{w) = w- VA-^ div w. 



In three dimensions, it is often convenient to rewrite equation (2.5) (by taking the curl 
of both sides) as 

d 

(2.7) — curl M + [m, curl m] = 0. 

The following was proved by the author in [?]. The main point is that the lin- 
earized Euler equation can be written in two ways, depending on whether one is right- 
translating or left-translating the perturbed tangent vector to the identity. For different 
purposes, either equation may be preferred. 



Proposition 2.1. The linearization of equation (2.5), obtained from a family u{t,a) 
of solutions by setting z{t) = da\^_f^u(t,a), is 

3z 

(2.8) — + P(V,n + Vuz) = 



with P given by (2.6) to keep z divergence- free. Writing z = rj^,v, we obtain the equation 
(2.9) p(^^^{Kv)+uj,xv^=Q, 

again with the consequence that v is divergence-free. Here A(t,x) = Dri(t,x)-''Dri{t,x) 
is the metric pullback and uJo{x) = curl'u(0,x) is the initial vorticity. 



In the study of fluid stability, the linearized Euler equation (2.8) is usually used alone 
to study Eulerian stability; however we can also study Lagrangian stability by incor- 
porating the linearization of the flow equation (2.4). From the geometric perspective, 
this corresponds to studying the Jacobi equation along the geodesic. The fact that we 
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can decouple the Jacobi equation this way was originally noticed by Rouchon [?]. It 
has been used to derive many properties of Jacobi fields and conjugate points; see for 
example [?]. 

Proposition 2.2. Ifr]{t) = exp^^{tUo) is a geodesic inV^{M) with velocity field u{t) = 
^ o rj{t)^^ , then the Jacobi equation along rj, right-translated to the identity, may be 
written as the pair of equations (2.8) and the linearized flow equation 

(2.10) ^ + V„i/- V,M = ^. 

Here divy = follows from the fact that div z = 0. 

The differential of the exponential map may be written as 

d{eWid)tuo{tZo) = yit) o ri{t), 



where y{t) is the solution of equations (2.8) and (2.10) with initial conditions y{0) = 
and z{0) = Zo- 

Corollary 2.3. We may also write the differential of the exponential map as 

d{^W^d)tuo{vo) = Dri{t)w{t), 
where w(t) solves the differential equation 

with initial conditions w{0) = and w'{0) = Vo- 

Proof. We have y{t) = ri{t)^w{t), and the equation ^ + [u,y] = z is equivalent (by the 



dt 

definition of the Lie derivative) to §i{^{t)^^y{t)) = r]{t)~^z{t). Thus we have ^ 
where v satisfies (2.9). The corollary follows. □ 

We know that the exponential map is defined and C°° on some neighborhood of the 
identity in P*(M) for s > n/2 + 1. It is not defined in the weak topology, since we do 
not have an existence and uniqueness theorem for solutions of the Euler equation with 
initial velocity. Thus the differential of the exponential map will only be defined 
only on some Vt C T^^V^^^M) Texp(fu<,)I^* (M). However, if Uo is sufficiently smooth, 
we can uniquely extend the differential to a continuous map on TidI^°(M) in the 
topology. (We emphasize that despite the notation, D°(M) is not a topological group 
or a smooth manifold in any known sense.) 

Definition 2.4. Let us denote by T^T'°(M) the closure of TrP^{M) in the topology. 
For any fixed a G and Uo G Tidl'^(M), we set E{a) = (dexpijj)^^^, and we define 
E{a): TidP°(M) Texp(a«„)2^°(M) to be the extension of E{a) to the L^-closure of 
Tidl^^(M). (This is a closed subspace of the space of all vector fields, by [?].) 

We can use E to extend the notion of conjugate points to the weak topology defined 
by the Riemannian metric. We say 77(a) is weakly monoconjugate to the identity if 
E{a) is not injective, and that 77(a) is weakly epiconjugate to the identity if E{a) is not 
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surjective. Similarly, 77(a) is strongly monoconjugate if E{a) is not injective, and ri{a) 
is strongly epiconjugate if E{a) is not surjective. 

The relations between weakly conjugate points are the same as those for conjugate 
points in strong metrics proved by Biliotti, Exel, Piccione, and Tausk [?], since those 
authors' proof relies only on the structure of the Jacobi equation rather than on global 
properties of infinite-dimensional manifolds. Thus we have: 

Theorem 2.5 (Biliotti et al.). We have the following for weakly conjugate points: 

• The set of weakly monoconjugate points is countable and dense in the set of 
weakly epiconjugate points; 

• Every weakly monoconjugate point is also weakly epiconjugate; 

• Any point that is strictly weakly epiconjugate (weakly epiconjugate and not 
weakly monoconjugate ) is a limit point of weakly monoconjugate points. 

These relations are not in general known to hold for strongly conjugate points in the 
present situation, where the Riemannian metric does not generate the if* topology 
in which we are working. However, we do have the following general results due to 
Grossmann [?]. These results rely only on the symmetry of the conjugacy relation. 

Theorem 2.6 (Grossmann). Weakly and strongly conjugate points both satisfy the 
following relationships: 

• Every (weakly /strongly) monoconjugate point is also a (weakly /strongly) epi- 
conjugate point. 

• // 77(a) is (weakly /strongly ) epiconjugate to the identity, and if the range of 
{E{a), E{a)) is closed, thenr]{a) is also (weakly /strongly) monoconjugate to the 
identity. 

In addition, the following relationships between weakly and strongly conjugate points 
are obvious: every strongly monoconjugate point is a weakly monoconjugate point (the 
kernel of the H'^ map is a subset of the kernel of the map), while every strongly 
epiconjugate point is a weakly epiconjugate point (the image of the map is a subset 
of the image of the map). The following proposition shows that for smooth geodesies 
in V^{M), strictly weakly epiconjugate implies strongly epiconjugate. For the basic 
facts on Fredholm operators that we will use, we refer to Taylor's Appendix A [?]. 

Proposition 2.7. Suppose the initial velocity Uo G Tif{DJ^{M) is C°°, so that the geo- 
desic T]{t) = expj^(tMo) is C°° as long as it exists. If ri{a) is strictly weakly epiconjugate 
to id, then ri{a) is strongly epiconjugate to the identity. 

Proof Assume, to get a contradiction, that 77(a) is strictly weakly epiconjugate but 
not strongly epiconjugate to id. 



We use the notation of Definition 2.4 By Proposition 2.6, 77(a) is strictly weakly 
epiconjugate if and only if the range of E{a) is not closed, which implies E{a) is not 
Fredholm. In addition, since 77(a) is not strongly epiconjugate to id, we know E{a) is 
surjective, which means it must also be injective; therefore E{a) is Fredholm. 
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Thus if we left-translate the Jacobi operators, as 

F{a) = dL.^a)-^ o E{a) : TidP^(M) ^ TmP^(M) 
F{a) = dLr,(a)-^ o E{a) : T,J)l{M) ^ T.^VUM), 

then -F(a) is also Fredholm, while F{a) is not Fredholm. If we consider in addition the 
operator 

where A = {1 + curl^)^/^, then since A'^ is an isomorphism between TidI'*(M) and 

TidI^°(M), we know F{a) is also Fredholm. The operators F{a) and -F(a) are both 
defined in the sa me sp ace, and thus we can consider their difference. 

By Proposition 2.3, the operators F{t) and F{t) both satisfy the differential equation 



(2.12) 



P 



,,,dF\ dF 



0, 



F(0) = 0, F'(0) = id. 



Thus the operator F{t) satisfies 

p( 



dt 



where A = (l+curl^)^/^. By considering commutators (and using the fact that [A, P] 
0), we have 



(2.13) P 



d f ^ , .dF\ 



-P 



Now since A(t) and Uo are smooth operators (for any t), we know that [A'^, A(t)] and 
[yl'^,ci;o] are both differential operators of order less than s, by the product rule; hence 
[A'^, A(t)]74~* and [A'', tUojA"'* are both compact operators on L^. Now given that F is 
a continuous operator from Ti^P^ C to itself, and given that P is a smooth operator 
from to Ti^V^ d can write 



(2.14) 



d 
dt 



d{F - F) 
Jt 



OJr. 



d{F - F) 
Jt 



K(t) 



where K{t) is a fixed compact operator from T\^T>^^ C to itself. Now F{Q) = and 

F(0) = 0; in addition -F'(O) = id and F (0) = id, so that F — F vanishes to both first 
and second order. 

Now consider the equation 

P j^{K{t)q{t))+u,q{t) =0. 

Since the operators are all bounded in L^, this is a linear differential equation in 
TidI^°(M), with solutions existing for all time; hence it has a bounded solution operator 
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G{t,T): TidV^iM) TidP°(M) such that q{t) = G(t,r)g(r) for all t, r G M. By 
Duhamel's principle, we can then write 

F{t) - F{t) = [ [ G{t, a) o K{(x) da dr. 
Jo Jo 

Since K{a) is compact for all a and G{t, a) is continuous, the composition is compact. 
The integrals are limits of sums of compact operators, hence also compact. 

Now by assumption we know that F{t) is Fredholm while F{t) is not Fredholm; this 
is a contradiction since the sum of a Fredholm operator and a compact operator must 
be Fredholm. □ 

We could also try to relate weak monoconjugacy to strong monoconjugacy, and this 
could be done most simply if we knew that every solution to (2.11) with w{0) = 



and w{T) = were actually as long as the flow and velocity field were C°°. (This 
would imply that any monoconjugate point corresponds to a zero of a C°° Jacobi 
field, and in particular to a zero of an H"^ Jacobi field.) While this is true in many 
cases, it is not always true, due to the fact that the order of a monoconjugate point 
may be infinite. 

The example of the 3-sphere from [?] is illustrative. In this example, we take a 
velocity field that is left-invariant under the standard group action on S'^; then it IS 
automatically a steady solution of the Euler equation and so generates a geodesic of 
1^fi{S^). Monoconjugate points occur at t = where m and n are any positive 
integers with m > n. Every such point has infinite order, and the infinite-dimensional 
space of Jacobi fields vanishing at both times is spanned by curl eigenfields, which are 
However, it is easy to find an infinite sum of such fields which converges in but 
not in for any s > 0, just by choosing the coefficients correctly. Hence there is a 
Jacobi field vanishing at t = and t = which is in but not in H'^ for any s > 0. 

In spite of this, it is still true that every weak monoconjugate point along this 
particular geodesic in Vfj_{S^) is actually a strong monoconjugate point, and we believe 
it is likely that this is always true. However the example above shows the proof may 
be somewhat subtle. 

The following theorem from [?] tells us when the first weakly conjugate point occurs, 
in terms of an ordinary differential equation. 

Theorem 2.8. Suppose t]: [0,T) — T>^{M) is a geodesic. Let u he the velocity field 
defined by ^{t,x) = u(t,ri{t,x)) . Let us define A{t,x) = Dri{t,x)'''Dri{t,x) and 
i-Oo{x) = curl'Uo(a^)- 

For each x G M, let r(x) > be the first time such that 

d / , , du\ , ^ du 
(2.15) -(a«.x)-)+.„(x)x-=0 

has a solution vanishing at t = and t = t{x). Then the first weakly conjugate 
point to the identity along rj occurs at mfr^izM t{x) . This point is either strictly weakly 
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epiconjugate (i.e., weakly epiconjugate but not weakly monoconjugate) , or it is weakly 
monoconjugate of infinite order. 



3. Locating epiconjugate points 
Since the Jacobi equation is closely related to the linearized Euler equation by Propo- 



sition 



2.2, we can use techniques developed by Lifschitz and Hameiri [?] and Friedlander 
and Vishik [?] for the latter. These authors showed that solutions of the linearized Eu- 
ler equation for sharply peaked initial data could be approximated near the peak by a 
certain ordinary differential equation, using a WKB approximation. Lifschitz-Hameiri 
proved a weaker estimate that is valid more generally than that of Friedlander- Vishik, 
but one that is sufficient for our purposes. (Friedlander- Vishik's estimate is valid only 
for steady solutions of the Euler equation on the flat torus T^, while Lifschitz-Hameiri 
is valid for all solutions in any geometry.) 

The basic technique of WKB analysis is to posit a solution in the form 

z = e'*/'^a + r, 

with e a small parameter, and expand in powers of e to obtain simple equations for $ 
and a. Then we prove that the error term r is bounded by e. The following computation 
and error estimate was performed by Lifschitz and Hameiri [?]. 

Theorem 3.1 (Lifschitz-Hameiri). Suppose u is a (possibly time- dependent) solution 



of the Euler equation (2.5) on a compact manifold M, on some time interval [0,T]. 
If z = e**/^a -|- r solves the linearized Euler equation 

dz 

(3.16) — + V„z + V^M = 2VA-Miv(V^u), div2 = 0, 

with initial condition z{0) = e**°/^ao, then for e small, the dominant terms satisfy 

(3.17) 9^+"^'^) = °' 
and 

(3.18) ^ + V,a + = 2gV<l', 

where 

(319) {JaU^ 

while the remainder term r satisfies 

(3.20) ||r(t)|U2 <C£ 



for all t G [0,T], where C is a constant that depends on T and derivatives of ^ and a. 
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Equations (3.17) and (3.18) look like partial differential equations, but they are more 
properly thought of as ordinary differential equations, which makes analyzing their 
solutions simpler. We can write them using right-translation as ODEs along a particular 
Lagrangian path, or alternatively using left-translation as ODEs in a particular tangent 
space. 



Proposition 3.2. Ifu is a solution of the Euler equation (2.5) with flow rj, then we can 
write the solution of equation ( |3.17[ ) as $(t) = ^o°vi^y^ withV^{t) = {Dr]{t)~^)*V(^o- 
In addition, we have for each x the formula a(t,ri{t, x)) = a{t) where a{t) is a vector 
field along t ^ r]{t, x) satisfying the equation 
(3.21) 



Da 



+ V„(t)M(t) = 2 



(V,(,)n(t),V$(t)) 



(v$(t),v$(t)) «(o) = «(0'-)' Ht),vm)^o. 

Alternatively, for each Xo we can write a(t,r]{t, x)) = Dr]{t, x)(^P{t)) , where P{t) G 
Tx^M for all t and P(t) satisfies 

(3.22) ^ {'K^^A{t)m) + {U ^0)^0 X f3{t) = 0, m = a(0, x), {m,Q = 0, 
where C,o = V$o(a;) and Uo = curlMo(x). 

Proof. The equation ^ + «($) = can be rewritten as ^ ($(t) or7(t)) = 0, which 
implies the solution $(t) = $0 o ri[t)^^ . The formula for V$(t) is a consequence of the 
chain rule. Finally the formula 



D 

dt 



da 



a{t,r]{t,x)) = — {t,r]{t,x)) + V„(t,,,(t,^))a(t, r7(t, x)) 



allows us to write (3.18) as the ordinary differential equation (3.21). 

To obtain (3.22) from (3.21), we apply (Dri)^ to both sides of (3.21) and obtain 



(3.23) 



Da 
~dt 



for some function c. The left side was computed in [?] to be, with a = Dri{j3), 



Da 
~dt 



+ Vq,m 



d 
dt 



(A(t)/3(t))+a;oX/5(t), 



using conservation of vorticity. Applying 7r£± to both sides of (3.23) gives (3.22). □ 



Corollary 3.3. Ifw{t) is a solution of equation (2.11) with initial condition w{0,x) = 
and Wt{0, x) = e**°(^^/^fo(x), then w{t, x) = e'*"*^"^^/^ /5(r, x) dT + f{t,x), where P{t,x) 
satisfies (3.22) with initial condition l3{0,x) = Vo{x) and ||f(t)||i2 < Ce fort G [0,T]. 
Here C a constant depending on T and derivatives of rj and Vo- 

Proof. We have tu(t, x) = JqV{t, x) dr and 



v{t) = v{t):'z{t) 



rj{t);'a{t) + rj{t);'r{t), 
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where a satisfies (3.18). Thus 

\\w{t) 



m\\ 



L2 



< 



< Ce 



< Ce 



Mr)-MT)\\dT 
t 

'Urr'Wc^ dr 



since rj is volume-preserving. With C = C /q ||^7(r) ^Hc-i dr, we are done. 
The main theorem of the paper follows. 



□ 



Theorem 3.4. Let rj he a smooth geodesic in P^(M) with Eulerian velocity field u. If 
for some Xo E M and some unit-length G T^^M , the equation 



(3.24) 



d 
di 



d'y 
~dt 



0, (7(t),Q = 



has a solution with 7(0) = and 7(a) = 0, then rj{a) is weakly epiconjugate to the 
identity along rj. 

Proof. Let 5 > be any number small enough that Riemannian normal coordinates 
exist in a ^-neighborhood of Xq- Choose normal coordinates {xi,X2,X3) with axes 
aligned so that = S^J,. and 7'(0) = d,j:.,\ , where 7 is the solution of (|3.24l) with 
7(0) = 7(a) = 0. Define Bs{xo) ^ M by ^o{xi,X2,X3) = Xi, and set = V$o = 
dxi + 0(5^). (We could extend $0 globally with a bump function, but this is not 
necessary, since we will only be working in this neighborhood.) Let M ^ M be a 
bump function with ip = 1 in [— |, = on M \ (—1, 1), and \ip'\ < 3 everywhere. 

Define a function a: M ^ M by a{xi,X2,X3) = X24' ^ ^^^^^^ + (and zero outside 

the coordinate neighborhood). Finally, define f o = ^ x Va. Then in the support of Vo, 
we must have Xi < 6, X2 < 5^, and 0:3 < 6. As a result, we have 



V7 2X1X2 , 



X2 I ^2 
1 ~'~ 3 



52 



+ 



+ 



Xr 



6^ 

x\ + X3 



2x3X2 
52 



X2 I ^2 
1 ~'~ 3 



+ 



Xr 



X2 
2 



2xo , / 



52 54 

and since X1X2 = 0{5^) and X3X2 = 0(5^), we have 



52 54 

,-yj2 I ^2 



a 



3:3 



/y>2 

X2 
54 



Va = 

Therefore, 

Vo(Xi,X2,X3 



/Y.2 I ^2 ^2 
"''I ~r -^a _|_ -^2 

52 54 



2x2 



54 



,-yj2 I ,y>2 ,y>2 

52 + 5? 



d,.,+0{5). 



/y«2 I ,-y.2 ^2 

"''I ~r -^3 _l_ -^2 
52 + 5^ 



^2 I ^2 ^2 
•^1 ~r -^3 _|_ -H 

52 ^ 54 



5,3+0(5). 
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Since X2 = 0(6"^), ip = 0(1), and ip' = 0(1), both terms in Vo are 0(1) in the support 
of Vo- 

Let us write T{t, x) : ^o{x)-^ C T^M — > ^o{x)-^ for the solution operator of 

d f ./ N d'-y , A (uJoix),Coix)) ^ , . d'l , 

- (.,,.,. A(t,.) H- LAJ^^„(,) X ^(M) = 

with initial conditions r(0,x) = and (9tr(0,x) = id, i.e., r(t,x)(xo) = 7(^) where 
7(0) = and 7(0) = Xo- Then T{t,x) is smooth and does not depend on 5; thus we 
can write r(a,x) = V{a^Xo) + 0(5). By assumption and our coordinate construction, 
we have V[a,Xo){dz) = 0. Since T(t,Xo) is an operator only in the space perpendicular 
to C,o = <9x.Jxo) and since r(a, Xo){drc:j) = by construction of our coordinates, we have 
for any x e L{xo)^ that T{a,Xo)ix) = {x, da:^)'^{a, Xo){d^^) + 0{5). 
In particular for the Vo constructed above, we have 

T{a,x){vo{x)) = {vo{x),d^,)T{a,Xo){d^,) + 0{6)=0{6), 

so that 

\\r{a,x){vo{x))\\L2=0{6)-Yol[snpp{vo)]=0{6'). 
On the other hand, we have ||^^o||l2 = 0(5"^) since Vo = 0(1) on supp(fo). 



By Corollary 3^, we have 

\\w{a)\\L2 < B5^ + Ce 

for some constant B. Since C depends on derivatives of f o, it may have a (5-dependence; 
however we are still free to choose £, and thus we can choose it small enough so that 
\\w{a)\\L2 = 0{6^). On the other hand, we still have ||w7'(0)||l2 = 0{6^). 

Now we consider a normalized sequence Wn of solutions with 5^ = ^ and e„ chosen 
so that 1 1 w„ (a) 1 1 2,2 = O(^) and (0)11^2 = 1. Then the series Yl'^=i'^n{o) converges 
to some divergence-free vector field in L^. We now have two possibilities: either E{a) 
is injective, or it is not. If E{a) is not injective, then T]{a) is weakly monoconjugate 
to id and thus also weakly epiconjugate to id, and we are done. Otherwise, if E{a) is 
injective, then we have Yl'^=i ^{^){'^'ni^)) convergent in L^, while Yl'^=i'^ni^) cannot 
converge in L^, so that E{a) does not have closed range. Hence ?y(a) must be weakly 
epiconjugate to id. Thus we are done. □ 



Remark 3.5. In the proof of Theorem 3.4, we assumed for simplicity that the dimension 
of M is 3. However the theorem holds generally as long as dimM > 2. In the more 
general case we need to replace the operator x ^ curlwo x x with x ^ ('•x'^'^o)^ where 
the generalized vorticity cross product is defined in arbitrary dimensions to satisfy (for 
any C) 

{{i^du'JX) = dxtixX) = (V^UoX) - {VcUo,x)- 
With this modification, all other parts of the proof are valid, with the expression xf + x"^ 

replaced whenever it appears with xf (in dimension two) or xl + x^ + xl + h (in 

dimension n > 3). 
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However, in dimension two the theorem happens not to give any useful information. 
This is because we must have (7,^0) = 0, which in two dimensions imphes that 7 = 
fC,^ for some function /. (Here ^j" is the standard rotation by 90°.) As a resuh, 
(iX^(7,7) = ffdXl{^^,^^) = regardless of the function /, so that equation (3.24) 
becomes 

|((e,A«.^-.)a|)-o. 

This equation clearly has no solutions f(t) vanishing at two distinct times, due to 



positive-definiteness of A, and hence Theorem 3^ does not yield any epi conjugate 
points at all in dimension two. 

This is as expected, since the exponential map is Fredholm in dimension two [?]. If 



there were nontrivial epiconjugate points obtainable by Theorem 3.4, it would imply 
(as shown in the next section) that epiconjugate points could occur in intervals along a 
geodesic; but Fredholmness of the exponential map implies that conjugate points must 
be isolated along a geodesic. 

On the other hand, the fact that the same technique works and yields genuine epi- 
conjugate points in any dimension higher than three implies that Fredholmness must 
also fail in any dimension higher than three. 



Although Theorem 3.4| allows us to locate many of the conjugate points along a 



particular geodesic, it will typically not get all of them. Essentially, the reason for 
this is the following: the first conjugate point is obtained from the solution of a three- 



dimensional second-order differential equation at some point, by Theorem 2.8, while 
the epiconjugate points detectable by the WKB method all come from the solution of 
a two-dimensional second-order differential equation. In essence, we are looking at the 
vanishing of the index form 

(3.25) /„„(-,,-,) = £(A(t,.)^.^l) + X -,(*). 5) A 

among vectors 7(t) vanishing at t = and t = t. Both problems have the same index 
form, but the WKB problem forces the vectors 'y{t) to all be orthogonal to the same 
fixed vector. Thus there are fewer solutions. 

Explicitly, suppose the first conjugate point location r along a geodesic is found from 



equation (2.15), such that that equation has only one solution vanishing at both t = 
and t = T, with nontrivial variation in all three directions. Then no two-dimensional 
vector solution can make the index form vanish, and thus the first epiconjugate point 



coming from Theorem 3.4 must be beyond the first conjugate point coming from The- 



orem |2.8[ We construct such an example in the following. 

Example 3.6. Consider the annular solid torus 5*^ x 5*^ x [a,b] (the region enclosed 
by two concentric tori) with coordinate system (x, y, z) in which the metric is given by 

ds^ = z dx^ + {dy + f{x) dz^ + dz"^, 
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for some periodic function f{x). If u = in this coordinate system, then curlw = dy, 
so that we automatically have [u, curl-u] = 0, which implies by (2.7) that m is a steady 
solution of the Euler equation. 

In the orthonormal basis ei = dx, 62 = dy, 63 = z^/'^{dz — f{x) dy), the operators 
in (2.15) take the form 



10 

A{t,x,y,z)= (0 1 k{t,x,z) 

k{t,x,z) l + k{t,x,zY 



and 



uJo{x,y,z) 




where k{t,x,z) = -/z(^f{t + x) — f{x)). The solution operator T of (2.15) satisfying 

T(0) = and V{0) = id is 

(3.26) 

(sint F(t) sint — G(t) cost cost — T 

-Git) t + J^k{T)'dT + J^[G{r)F'{T)-F{T)G'{T)]dr -Fit) 
1— cost —F(t) cost — G(t) sin t sint 

where F(t) = k^r) cost dr and G{t) = k^r) sinr dr (suppressing the dependence 
on X and z for simplicity). 

There is a conjugate point at t iff there is some vector Vq with T(t)(fo) = 0, i.e., iff 
det T(t) = 0. We can easily compute 



detT(t) = - sint F(t)2 + G(t)2 
+ 2(1 - cost) 



t+ [ k{r)UT+ [ [G{T)F'{T)-F{r)G'{r)]dT 
Jo Jo 



For a particular example, e.g., f{x) = sinx, it is easy to check numerically that for any 
(x, y, z), the first time where det T(t) = is t = 27i. (This is probably true in general, 
but it is not important here). 

The unique (up to rescaling) vector Vo for which T(27r)(w„) = is 




and the corresponding solution to (|2.15|) is 

l{t) = 



F{2ti) sint + G'(27r)(l - cost) 

G{2Tx)F{t) - F{2n)G{t) 
F(27r)(l - cost) - G(27r) sint. 



There is no vector C,o for which {'-/{t),^^) = for all t. As a result, the reduced 
index form does not vanish at t = 27r, and so the WKB method does not predict the 
epiconjugate point at t = 27r. 



THE WKB METHOD FOR CONJUGATE POINTS 



15 



4. Conjugate point intervals 



We now explore the consequences of Theorem 3.4 First, we have an easy theorem 
on continuity of conjugate point locations with respect to the fixed initial vector C,o- 

Theorem 4.1. If for some Xo & M and some unit vector ^ T^^M , the equation 
(3.24) has a solution -yit) with 7(0) = and •yito) = 0, then for any sufficiently 
close to C,o, there is a solution ^(t) with 7(0) = and 7(^0) = for some to- 
Proof. Fix an orthonormal basis of vector fields ei, 62, and 63 in a neighborhood of Xo- 
Using spherical coordinates 6 and 0, we can write for any x near Xo 

= sin 9 cos 0ei | ^ + sin 9 sin 0e2 1 ^ + cos 9e3 1 ^ 

^1 = cos 9 cos 0ei I ^ + cos 9 sin 0e2 1 _^ — sin 6^63 1 ^ 

^2 = -sin0ei|^ + cos0e2|^. 

Then any vector orthogonal to C,o must be a linear combination of ^1 and ^2, so that 
the equation (3.24) becomes, with u{t) = + g(t)^2, 



d 

^'■''^ dt 



(^i,A(t,x)ei) (ei,A(t,x)6)\ ff{t) 
(^i,A(t,a;)6) {^2,A{t,xK2)J \m 



+ Mx) X ^1,^2) 



A (m\ ^ (0 

-1 0) {git)J VO 



In this way, we can consider the dependence of the conjugate time to on 9, 0, and x. 
By the general theory of oscillation for self-adjoint systems (see Reid [?]), equation 

(4.27) has a solution satisfying m(0) = and u{to) = if and only if the index form 

(4.28) I^(f^,+g^,J^,+g^,) = 

(An(t, 9, <P)f{tf + 2Ai2(t, 9, ^)f{t)g{t) + A22(t, 9, 4>)g{tf + 2uo{9, <P)f{t)m) dt, 

defined for vector functions u{t) vanishing at t = and t = T, is negative on some 
subspace for T > to and positive-definite for T < to- (Here, for each x G M, we have 
kj{t,9,<P) = {^i,A{t)Q and Uo{9,<P) = K x ^1,6)-) 
In fact, if we set 

Jo llitWdt 

then for any fixed 9o and (po, we will have for any T < to that the infimum sat- 
isfies inf^ 7(6^0, 0o) (7) = X{T,9o,(f>o) > 0, while if T > to, then inf^ 7(6^0, 0o) (7) = 
A(T, 9o, (po) < 0. We can prove this by noting that, by standard Sturm-Liouville theory 
for self-adjoint systems, J [9, 4>){'y) attains a minimum at a certain eigenvector field 7(t) 
of a self-adjoint operator, and the minimum is an eigenvalue which we can denote by 
X(T,9,(f)). This eigenvalue depends continuously on T, 9, and (p by the usual theory. 
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As a result, for any fixed T < to, if 6 and are sufficiently close to 6o and 0o, we will 
have A(T, 6,(f)) > 0; and for any fixed T > to we will have A(T, ^, 0) < 0. As a result, 
by continuity with respect to T we must have X{to, 6,(j)) = for some to- □ 

The theorem above gives existence of (possibly multivalued) functions t of spherical 
variables. For any fixed 6o and 0o we have an open domain in the sphere, containing 
9o and (po, on which t{9,(j)) is defined. Hence t{9,(f)) is defined on some open subset 
of the 2-sphere. We can understand its behavior better by examining the differential 
equation directly, as follows. As a consequence, we will prove that t is a differentiable 
function of 6 and 0. 

Proposition 4.2. At a point Xo & M and for a unit C,o £ T^^M with {C,o,^o{xo)) 7^ 0, 
the equation (3.24) has a solution with 7(0) = andj{a) = if and only if the solution 
of the 2x2 matrix equation 

dW 

(4.29) 



dt 



+ {uJo{Xo),Qj{7T^±A{t,Xo)n^±) 



with W{Q) = id satisfies TrW{a) = 2. Here J 



0-1 
,10. 



Proof. Let us write the solution 7(t) of (3.24) with 7(0) = and 7'(0) = f o as a matrix 
operator u{t) = S(t){vo), with S(t) : Then the 2x2 matrix T must satisfy 



(4.30) 



d 



dt 



dS\ 



dt J 



dS 
lit 



0, 



with ^(0) =0 and f (0) = id. 

Writing W = tt^±A ^ and using the fact that A(0) = id, we have the initial condition 
W{0) = id, while W satisfies (4.29). In addition, since vr^Avr^ is symmetric and J is 
antisymmetric, we have Tr ^W~^ = 0, so that det W{t) = 1 for all t. On the other 



hand, if we integrate (4.30) in time, we obtain 

W{t) + {UoOo)JS{t) 



id. 



Since by assumption {C,o,^^o) 7^ 0, there is a f o G with S{a){vo) = if and only 
if W{a) has 1 as an eigenvalue. Since is a 2 x 2 matrix with det W = 1, the 
eigenvalues of W satisfy — Tr W{a)\ + 1 = 0, so that A = 1 is a solution if and only 
if Tr W{a) = 2. □ 

Now to study differentiability, we fix the dependence of the solution on G S"^ as 
in Theorem 



4.1 



(4.31) 



then equation (4.30) becomes 
dW{t,e,(f)] 



dt 



jQ{t,e,(f))w{t,e,(i)) = 0, 



with 9(t, 9, 0) = {uJo{xo),^o){7^^^A{t, Xo)'n'^±) ^. In this way, we get a function W{t, ^o) 
defined on the 2-sphere. 
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Theorem 4.3. Let W{t,^o) be the solution of equation (4.31) with W{0,^o) = id, and 
suppose we have TrW(to,^o) = 2 for some to > and ^ S"^ with {^o,^o{xo)) 7^ 0. 
Then we can solve the equation TrW{t,C,o) = 2 for t in terms of ^, differentiably in 
any neighborhood of C,o- 

Proof. Let us consider the 2-sphere as a manifold; it is enough to prove that if ?/ is a 
coordinate on S'^ in some system, then ^ exists. So let us suppose the other coordinate 
as fixed, and think of W{t,y) as depending only on y. We assume y = yo corresponds 
to the specified point in S"^- 

Since det W{t,y) = 1, we can show using the Cayley-Hamilton Theorem that we 

have [W{to,yo) - id] = 0. Thus W{to,y) - id = for some nilpotent matrix A^, 
which must look like = cJC,C^ for some unit vector and real number c. Directly 
from equation ( 4.31[ ), we have 



^ Tr Wito, yo) = - Tr ( je(t„ yo) (id + cJC(^)) 
= -cTT{jQ{to,yo)JCe) 

= -c((c, Je(t„,y„)JCC^(C)) + {JCJe{to,yo)JCe{JO) 

= c{JC,Q{to,yo)JO- 

In the first line, we used the familiar fact that if A is antisymmetric and B is symmetric, 
then Tt AB = 0. (We will use this again many times throughout the proof.) Now the 
term {J(,Q{to,yo)JO is positive since 9 is positive-definite, so that if c 7^ 0, then 
'^'^'^^l*"'^"-' 7^ 0. As a result, we can solve for t in terms of y by the implicit function 
theorem (and hence obtain j^) if W{to,yo) 7^ id. On the other hand, if W{to, yo) = id, 
then ^^L^ {to, yo) = 0. Thus we have to work a bit harder. 



So now we suppose that W{to, yo) = id. Since the coefficients of (4.31 ) are smooth in 
y, we can differentiate W{t,y) with respect to y. If we differentiate (4.31) with respect 
to y and write ^ = WM, then M must satisfy the equation 

and since Tr{W-^J^W) = Tr(jg) = 0, we have | Tr M(t, y) = for aU t and y, so 
that TrM{t,y) = for all t and y. Thus since W{to,yo) = id, we have ^'^gj^ {to, yo) = 
TTM{to,yo)=0. 

So both partial derivatives of TrW are zero at {to,yo), and thus we can look at 
second derivatives. If ^ exists, then we must have (by differentiating Ti W {t{y),y) 
twice with respect to y) 

d^TiW, ,fdtV d'^liW , ,dt d'^TiW , , ^ 
-^{to,yo) ^^J + 2^^(^0,1/0) ^ + ^^(to, 1/0) = 0, 
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and in order to be able to solve this for ^, we must have 



(4.33) 



{to,yo) 7r7:^{to,yo) a 9 {to,yo) > 0. 



dtdy 



dy"^ 



Conversely if this condition is satisfied, we obtain either one or two values for ^, 
corresponding to a crossing of at most two solutions. 

Observe first that, quite generally, if is a matrix with det W = 1, then for any 
parameter r we have Tr(iy~^^) = 0. If we then differentiate with respect to any 
other parameter s, we have 



Now in our special case, we have W{to,yo) = id so that 



(4.34) 

while 

(4.35) 

and 

(4.36) 



dOdt 



dy"^ 



(to, yo) = Tr 



/dW 

{to,yo) = Tr -7—{to,yo) 



[to,yo) -7^{to,yo, 



dt 



(to, yo) 



V dy 



V dt 



dy 



fdW 

-2 det y-^{to,yo] 



{to, yo/ 



-2 det 



dy 
(dW 



[^-Q^{to,yo) 



(these last two are consequences of the Cayley-Hamilton theorem) 

to 

'dW dW 



Thus the inequality (4.33) translates into 

2 



Tr 



dW dW 



dt dy 



4 det 



dt dy 



> 



at {to,yo)- 



Now the equation Tr — 4 det A > for a 2 x 2 matrix A is equivalent to the require- 
ment that A has real eigenvalues. Thus we will have (4.33) if 



(4.37) 



dW 

~dt 



dW 



dy 



dQ 



ito,yo)^ito,yo) = JQito,yo) W{T,yo) ' J —{T,yo)W{r,yo) dr 







dy 



has real eigenvalues. To establish this, we eliminate the J matrices using the easily- 
proved formula JZ = (det Z){Z'^)~^J along with = —id to obtain 



dW dW 

-{to, yo) -^{to, yo) = det (6 (to, yo))Q{to, yo) 



to QQ 

W{T,yo)^-^{to,yo)W{T,yo) dr. 



r\, \"0, ao; rs 

dt dy 

We now have the desired matrix as a product of a symmetric positive-definite matrix 
and a symmetric matrix. It is easy to see that such a product must have real eigenval- 
ues: let C be a symmetric positive-definite matrix and D a symmetric matrix. In an 
eigenvector-basis of C, we have C = ( o ° ) and D = ( g ^ ) , so that CD = ( ) • The 
expression [Tr (CD)]^ — 4det (CD) becomes (Aa — yUc)^ + 4Ayu6^, and this is nonnegative 
since \fi > 0. So CD must have real eigenvalues, and this establishes ( |4.33[ ). □ 
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Corollary 4.4. For each point Xq G M, there is a family of functions /(^o) defined 
on open subsets of the 2-sphere, such that for every there is a weakly epiconjugate 
point to the identity located at t = f{C,o)- As a result, the set of all weakly epiconjugate 



point locations obtainable by the technique of Theorem \3.4\ consists of the closure of the 
union of open intervals. 

Proof. For every component of the domain of /, the image of /. So we have a family 
of intervals (possibly open, closed, or half-open) as ranges over S'^. The actual set of 
epiconjugate points obtained this way must be closed, since if we have an sequence of 
conjugate point locations then any limit point must be an epiconjugate point location. 
Hence we obtain some union of closed intervals. □ 



By Theorem 2.5| and Theorem 2/7, we can conclude the following about the intervals 



obtained from Corollary 4.4 



Corollary 4.5. All of the points in the nontrivial intervals obtained by Corollary \4. 4 
are strongly epiconjugate. 

Proof. Take any nontrivial closed interval of weakly epiconjugate point locations [a, b] C 



M"*". By Theorem 2^, there is a countable dense subset of this interval consisting of 
weakly monoconjugate points. All other points in the interval must therefore be strictly 
weakly epiconjugate, and thus strongly epiconjugate. Finally, since a hmit of strongly 
epiconjugate points is also strongly epiconjugate, the entire interval must consist of 
strongly epiconjugate points. □ 

It is not clear whether this technique actually generates more than one interval; in 
the simple examples one can work out explicitly, there is always just a single interval 
extending to infinity. The following is a simple yet seemingly typical example. 

Proposition 4.6. Suppose u is a Killing field on M . (Such a u is a steady solution of 
the Euler equation, by Misiolek The epiconjugate points along the corresponding 
geodesic rj form an interval [r, oo) for some r > 0. 

Proof. Since m is a Killing field, the metric puUback A{t, x) must be the identity for all 



t and X. Thus the equation (3.24) takes the form 



with 7(t) orthogonal to for all t. It is easy to see that this equation has solutions 
vanishing at times 

{uJo[Xo),t,o) 

for any unit vector G Tx^M. The infimum t{xo) = inf^^g52 t{xo,C,o) occurs when C,o 
is parallel to ujoixo) and is t(xo) = i— The absolute minimum over all Xo G M is 

inf t(xo) 
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the time r(xo, ^o) 



On the other hand, as approaches a vector orthogonal to cOoixo 
obviously approaches infinity for any Xo- 

Now it is also not hard to see that the actual first conju gate point occurs at r 
(this is proved in [?]). As a result, in this case Theorem 



2lT 



point locations 



3.4 



predicts all conjugate 
□ 



The behavior seen above for Killing fields appears to be typical; for the fiows whose 
conjugate points one can work out explicitly, epiconjugate points occur in intervals 
extending to infinity. It is conceivable that for nonsteady fiows, the metric pullback 



A may increase fast enough with time to prevent solutions of (3.24) from vanishing at 



large times. Obviously, if the solution of the Euler equation blows up at a finite time, 
then all bets are off. 
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